Dynamical Casimir effect for magnons in a spinor Bose-Einstein condensate 
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Magnon excitation in a spinor Bose-Einstein condensate by a driven magnetic field is shown 
to have a close analogy with the dynamical Casimir effect. A time-dependent external magnetic 
field amplifies quantum fluctuations in the magnetic ground state of the condensate, leading to 
magnetization of the system. The magnetization occurs in a direction perpendicular to the magnetic 
field breaking the rotation symmetry. This phenomenon is numerically demonstrated and the excited 
quantum field is shown to be squeezed. 

PACS numbers: 03.75.Mn, 03.70.+k, 42.50.Lc, 42.50.Dv 



I. INTRODUCTION 

Vacuum fluctuations play an important role in a vari- 
ety of situations in quantum physics. For instance, the 
static Casimir effect [lj originates from vacuum fluctu- 
ations in the electromagnetic field and in the electronic 
states in matter. If the definition of the vacuum state 
depends on time due to a time-dependent external con- 
dition and the system cannot follow the instantaneous 
vacuum state adiabatically, the vacuum fluctuation ma- 
terializes as real particles @, d, 0, [BJ . This phenomenon 
is called the nonstationary or dynamical Casimir effect 
(DCE). 

The DCE has been extensively studied [y], especially 
for photons and the massless scalar field. When the mir- 
ror of an optical cavity is rapidly moved 0] or the di- 
electric constant of the matter in a cavity is rapidly al- 
tered @, 0, 03, EH, Photons are created in the cavity 
even if the initial state of the electromagnetic field is 
in the vacuum state. The generated photons are in a 
squeezed state 12J, which modifies the Casimir force ex- 
erted on the mirrors [l3| . The finite temperature correc- 
tion to the DCE and decoherence via the DCE [l5j 
have also been studied. Since there is no intrinsic Hamil- 
tonian in the original formulation of the moving mirror 
problem 3 , the effective Hamiltonian approach has been 
developed [M E3, El, El • 

We have derived an effective 
Hamiltonian for the moving mirror problem by quantiz- 
ing both the electromagnetic field and the polarization 
field in the mirrors [20]. 

However, photon creation by the DCE has not yet been 
experimentally observed. This is because the mirror of 
a cavity must be vibrated at a frequency of the order of 
GHz (at least for a microwave cavity) in order to reso- 
nantly amplify the photon field [2l[ . Recently, the INFN 
group proposed a method to detect the DCE using 
a semiconductor layer illuminated by laser pulses, which 
enables the rapid displacement of the position of the cav- 
ity mirror. The proposed experiment using this method, 
however, has not been completed to date. 

In the present paper, we propose a novel system to 
realize the DCE: a Bose Einstein condensate (BEC) of 
an ultracold atomic gas with spin degrees of freedom. 



In this system, magnons are created through the DCE. 
The time-dependent external condition that causes the 
DCE corresponds to the time-dependent magnetic field 
applied to the BEC. Unlike in the case of photons, the 
typical energy of magnons in a BEC for a magnetic field 
of ~ 1 G is ~ h x 100 Hz (h: Planck constant). Magnetic 
field modulation at this frequency is experimentally feasi- 
ble. Continuous amplification of magnons by an oscillat- 
ing magnetic field leads to magnetization of the system, 
which can be observed by in situ measurements [231 ] . 

In a broad sense, quasiparticle excitations in a nonsta- 
tionary BEC may be regarded as the DCE. For exam- 
ple, a time dependent trapping potential [24j . a rapid in- 
crease in the interatomic interaction 12511 , and collapse of 
a BEC by an attractive interaction [26[ generate Bogoli- 
ubov quasiparticlcs. In these cases, the BEC itself is also 
excited and its shape depends on time, whose dynam- 
ics is described by the mean-field Gross-Pitaevskii (GP) 
equation. However, it is difficult to distinguish the quasi- 
particle excitation from the mean-field excitation, and 
therefore these systems are unsuitable for demonstrating 
the DCE. In contrast, in our model, the time-dependent 
magnetic field excites only the vacuum fluctuation in the 
initial quasiparticle vacuum state, giving an ideal testing 
ground for the DCE. 

This paper is organized as follows. Section [TT] formu- 
lates the problem. Section IIIII discusses the relation of 
the magnon excitation in the present system with the 
DCE. Section HVl numerically demonstrates the proposed 
phenomena using the GP equation with quantum fluctu- 
ations. Section fVl provides discussion and conclusions. 



II. FORMULATION OF THE PROBLEM 

A. Hamiltonian for the system 

We consider spin-1 bosonic atoms with mass M con- 
fined in an optical trapping potential U(r). The single- 
particle part of the Hamiltonian without a magnetic field 
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is given by 



B. Bogoliubov approximation 



H = J2 J dr$l(r) 



4> m (r), (1) 



where ip m (f) is the field operator that annihilates an 
atom with spin magnetic quantum number m — —1,0,1 
at the position r. 

The interatomic interaction for ultracold spin-1 atoms 
is described by the s-wave scattering lengths ao and a 2 , 
where the subscripts and 2 indicate the total spin of 
two colliding atoms. The interaction Hamiltonian can 
be written in spin-independent and spin-dependent parts 
as [27L|28j 



Hint = / dr c : p 2 (r) : +ci : F(r) ■ F(r) 



(2) 

where the symbol :: denotes the normal ordering and 
1 

P(r) = £ 4>l(r)$ m (r), (3) 

rn — — 1 

^(r) = E $l(r)(f)m,m4m'(r), (4) 

m,m' 

with / = (f x , fy, f z ) being the vector of the spin-1 3x3 
matrices. The interaction coefficients cq and c± in Eq. @ 
are given by 



co 



AitTi ao + 2a 2 
~M 3 

a 2 — ao 
~~M 3 ' 



(5) 
(6) 



We restrict ourselves to the case of the hyperfine spin 
F = 1 of an alkali atom with nuclear spin 7 = 3/2 and 
electron spin 5=1/2 (e.g., 23 Na and 87 Rb). Because of 
the hyperfine coupling between the nuclear and electron 
spins and their different magnetic moments, the Zccman 
energy is a nonlinear function of B [29( . Taking the first 
and second order terms, the Hamiltonian becomes 



dr 



Pl B(t)F z + p 2 B 2 {t) + 



(7) 

where p\ and p 2 are the linear and quadratic Zeeman 
coefficients, respectively, and we assume that the uniform 
magnetic field B(t) is applied in the z direction. We 
define the quadratic Zeeman energy as 



q(t)=p 2 B 2 (t), 

which is positive for the F = 1 hyperfine state. 
Thus, the total Hamiltonian has the form 



H = H + H l 



flh 



(8) 



(9) 



Since the linear Zeeman term in Eq. ([7]) commutes with 
the other part of the Hamiltonian and only rotates the 
spin uniformly around the z axis, we neglect the linear 
Zeeman term henceforth. 



The initial state considered in the present paper is the 
ground state for an initial value of the quadratic Zeeman 
energy, q = g(0), under the restriction 



dr(F z (r)) = 0. 



(10) 



In Sec.|V]we discuss how to generate this state. The spin- 
independent interaction coefficient Co must be positive for 
the existence of the ground state. Either for c\ > or 
for ci < and q > 2|ci|p, where p is a typical atomic 
density, almost all atoms are in the m — state for this 
initial state [3(| ■ We therefore employ the Bogoliubov 
approximation by setting 



■0±i(r) 



^o(r) + Mr) 



(11) 
(12) 



where v I / o( r ) is a real function that minimizes the energy 
functional, 



-En 



dr 



2M 

fi is the chemical potential, 



dr 



2M 



V + u *o + c n 



(13) 



(14) 



and (f> m (r) are the fluctuation operators. The wave func- 
tion ^oif) is normalized as / dr^o( r ) — N with N being 
the number of atoms. The Heisenberg equations of mo- 
tion for 4>m are then written as [271 



dt 



-— V 2 + [/-/ij 0o+c o ^ 2 (20 o 



>±1 



>±1 



(15) 
(16) 



where we neglect the second and third orders of 
4> m . Equation (|15[) is identical with that of a scalar 
BEC, which gives the Bogoliubov spectrum E^ — 
\/ £ fc( e fc + 2co x I'o) f° r a uniform system, where = 
Ti 2 k 2 /{2M). 

If q is constant in time, we can solve Eq. (I16|) by setting 



J2 [u,(r)b ± , x e-^ lh ^ + v* x (r)bl x e iE ^ 



(17) 

where b±,\ are bosonic operators for quasiparticles and 
U\, v\, and E^' are determined by the Bogoliubov-de 
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Genne equations, 



2M 



2M 



X7 z + U + q + c *a + ci*o - M ) u A + ci*5ua 

(18) 

V 2 + U + q + C * 2 + Cl * 2 - /i ) V X + Ci* 2 M A 



(19) 



The quasiparticles created by b± x are excitations of m — 
±1 states, which we call magnons. 

For a uniform system, Eqs. (fT8| and (fT9|) are solved to 
give 



''fr 




(20) 
(21) 



where V is the volume of the system and 

E { h 1] = y/(e k +q)(e k + q + 2c 1 **). (22) 

When all Bogoliubov energies are real, the system is dy- 
namically stable. If ci < and q < 2|ci|\I>o, E( l- d22J) 
is imaginary for long wavelengths and the system is 
dynamically unstable against spontaneous magnetiza- 
tion [3lL [32I ]. In the present paper, we consider such 
parameters that the system is dynamically stable. 



C. Effective Hamiltonian for oscillating q 

Here we assume that the quadratic Zeeman energy q(t) 
oscillates as 



q(t) = <2o ( 1 + S q sin Qt) 



(23) 



through the oscillation of the strength of the applied mag- 
netic field. For 8 q = 0, the Bogoliubov Hamiltonian for 
magnons is given by 



(24) 



where E x , 6+,a, and S_.a are energies and annihilation 
operators of magnons defined at q = qo. The magnon 
Hamiltonian for S q 7^ is then written as 



H« = H?=l° + q Q 6 q sin Qt / dr ^ + ^-i 



When the last term in Eq. (|25|) is treated as the pertur- 
bation, its interaction representation becomes 



S q sin rit / dr 



(U* X U X > +V* x Vy) 



bl, x b-,X' 



-e-^ +Ei >^(u x v x ,+v x u x/ )b + . x b^ x/ 



-H.c. + const. 



(26) 



-,-fc i 1 



where H.c. denotes the Hermitian conjugate of the pre- 
ceding term. If the frequency Q is resonant with (e{ + 
E^yn for some A and A', the terms proportional to 
£>+,a&-.a' and V + yp_ y become dominant. 

When the system is uniform and f2 is resonant with 
2Ej^' /h, the Hamiltonian in Eq. (|2"6"|) reduces to 

H cS = -iqo5 q V^T \u k v k \ (b +tk b--k ~ b f +yk b[ 

k 

(27) 

where denotes that the summation is taken for 

h 2 k 2 /(2M) ~ 2E ( k 1] . The form of Eq. (J27J) suggests that 
the magnon state develops into the squeezed state by the 
oscillation of q, which will be confirmed in Sec. |IVB] nu- 
merically. 



III. RELATION WITH THE DYNAMICAL 
CASIMIR EFFECT 

We now discuss the relation of the present system to 
the usual DCE. 

First let us review the DCE for an electromagnetic 
field. The initial state is the vacuum state for some static 
configurations of mirrors and dielectrics. The initial ex- 
pectation values of the electric and magnetic fields there- 
fore vanish: (E) = (H) = 0. These operators obey the 
Heisenberg equations of motion, 



d_ 

Of 



e(r,t)E{r,t) = V x H(r,t), 



(28) 



d_ 

df 



H(r,t) 



-V x e(r,t)E(r,t), (29) 



(25) 



where e is a dielectric constant, and the electric field op- 
erator E must vanish at the mirrors. From Eqs. (|2"5]) 
and (|29[) . (E) and (if) always remain zero as in classi- 
cal electrodynamics, whereas the vacuum fluctuation of 
the electromagnetic field can be amplified, leading to the 
creation of photons. 

For the present spinor BEC system, the initial state 
is assumed to be the ground state satisfying Eq. (flQ)) . 
Since the Hamiltonian ^ has spin-rotation symmetry 
around the z axis, the expectation value of the transverse 
magnetization for the initial state vanishes: (-F+) = 0, 
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where F + — F x + iF y . In the Bogoliubov approxima- 
tion, using Eqs. and (HH), the magnetization opera- 
tor F + = v / 2(V'iV'o + ^W-i) reduces to 



F + 



(30) 



From Eq. ([16]), the Heisenberg equations of motion for 
re ol 

dF+ 



F + are obtained as 



dt 

dfl+ 
~df 



+ (c + 2ci)*2 _ p 



F+, 



where we define 



n 4 



n+j(3i) 



(32) 



(33) 



The expectation value of this operator also vanishes, 
(LT+) = 0, for the initial state. We find from Eqs. (pTT]) 
and (|32| that (F+) and (11+) always remain zero. On 
the other hand, their quantum fluctuations can be am- 
plified due to the temporal variation of q(t), leading to 
the creation of Bogoliubov quasiparticles, i.e., magnons. 
Thus, the present situation is similar to that of the DCE 
in the electromagnetic field in that the amplification of 
quantum fluctuations and (quasi)particle creations occur 
while the expectation values of the fields remain constant 

((F+) = (n+) = o). _ 

What kind of excitation phenomenon can we regard 
as the DCE? For example, can the excitation of rip- 
ples on water in a vibrating bucket be attributed to 
the DCE? The answer is partially yes, because ampli- 
fication of the quantum fluctuation on the water surface 
does occur, which creates ripplons. However, this quan- 
tum excitation is overwhelmed by the classical excitation 
of the water surface and identification of the quasipar- 
ticle excitations should therefore be extremely difficult. 
Thus, the suitable condition for studying the DCE is that 
only the quantum fluctuation is excited and the classi- 
cal field, i.e., the expectation value of the relevant quan- 
tum field, remains constant during the temporal vari- 
ation of external parameters. For the electromagnetic 
field, (E(r,t)) = (H(r,t)) = 0, even when the position 
of the mirror and e(r,t) of matter are changed. The 
present spinor BEC system is also suitable for studying 
the DCE, since (F+(r,t)) = (f[ + (r,t)) = holds during 
the change of q(t). 



IV. NUMERICAL ANALYSIS 

A. Mean-field theory with quantum fluctuations 

In this section, we numerically demonstrate the DCE 
of magnons in a spinor BEC. Since performing a full 



quantum many-body simulation is difficult, we employ a 
mean-field approximation taking into account the initial 
quantum fluctuations. Recently, this method was used to 
predict spin vortex formation through the Kibble-Zurek 
mechanism in magnetization of a spinor BEC [33j. 

The mean-field GP equations for a spin-1 BEC are 
given by 



dtp 



±i 



dt 



-—^ + V + q(t) + c oP 



^±1 



+ci ( -^Ft^o ± F z ip ±1 



dt 



in -ET = -7T77 V +V + c p)^ 



2M 



-^(i^i+i^V-x). 



(34) 



(35) 



where tp m (r,t) are the macroscopic wave functions and 
p, F z , and F± = F x ± iF y arc defined by the forms in 
Eqs. ([3]) and ((4]) in which tp m are replaced by ip m . The 
initial wave function for tpQ is the ground state solution 
of Eq. (fl~3|) , which is a stationary solution of Eq. |35|) 
for ip±i = 0. Although (F±) = for the exact many- 
body initial state, we do not set ip±i = 0, since the right- 
hand side of Eq. (|34| vanishes and no time evolution is 
obtained. 

We include a small initial noise in ij)±x to reproduce 
vacuum fluctuation in magnetization within the Bogoli- 
ubov approximation. To this end, we consider the mag- 
netic correlation function for the Bogoliubov vacuum, 

<F + (r)F_(r')> * 2* (r)*o(r') (r)tf(r'), (36) 

A 

where f\(r) = u\(r) + V\(r) and we have used Eqs. (fT7| 
and pop . To find the appropriate mean-field initial states 
of ip±i, we assume the form 



lM»0 = E Mr)&±,A + vl(r)b^ x ] , (37) 

A 



where b± \ are random numbers whose probability dis- 
tribution is determined below. Substituting Eq. (|37[) into 
F+(r)F-(r') and taking the average with respect to the 
probability distribution, we obtain 

(F+(r)F_(r'))av g - 2* (r)* (r') 

xE[/AM/A<(0<& + ,Ak + ,A'}av g 
A, A' 

+/A(r)/ A */(r / )(6-,A6!L,v) W8 
+f\(r)f x .{r'){b-,xb+, x >) mK 
+mr)mr')(b* + ,xb*_^U s ]-(38) 

Using the fact that both (u, v) and (u*,v*) are solutions 
of Eqs. (fT5)) and (fT§)) and assuming that the random vari- 
ables 6+^a and t>_ x obey the same distribution, we find 
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that Eq. (f38|) coincides with Eq. 
distribution satisfying 



36)) by a probability 



(6 



±.A°±,A' 



(b 



±,A°=F ,A'/avg 



= 0. 



(39) 



We also assume (&±,A&±,A')av g = (6±,A^,A')avg = 0. 
Thus, the initial mean-field states (|37|) with probabil- 
ity distributions obeying Eq. (|39p reproduce the same 
magnetic correlation as for the Bogoliubov approxima- 
tion. Similarly we can show that the time evolution of 
(F + (r)F_(r'))av g also agrees with that of (F+(r)F_ (r')> 
at the level of the Bogoliubov approximation. 



Numerical results 



We numerically solve Eqs. (]34[) and (|35|) for the time- 
dependent magnetic field that gives a sinusoidal oscil- 
lation of the quadratic Zeeman energy as in Eq. (|23p . 
For simplicity, we consider a uniform two-dimensional 
(2D) system with a periodic boundary condition. We as- 
sume that 10 7 23 Na atoms are in a 100 /im x 100 /im 
space with atomic density n — 2.8 x 10 14 cm~ 3 . In a 
recent experiment [35| . 02 — ao oc c\ for an F = 1 23 Na 
atom was precisely measured to be about 2.47 Bohr radii. 
The system is therefore dynamically stable according to 
Eq. I]22p. The initial state is ipo = \pn anc ^ V'ii are given 
by Eqs. (|37[) and (|39p in which the complex random vari- 
ables b± ! k are assumed to follow the Gaussian distribu- 
tion P(b) — y / 2/-Ke~ 2 ^ 2 . The initial random variables 
b±jf are cut off for 2ir/k < 1.5 /zm. The strength of the 
magnetic field at t = is 500 mG, which corresponds to 
q Q ~ ft, x 69 Hz ~ 1.02ciTi. 

Figure Q] (a) shows time evolution of the average 
squared transverse magnetization [34j , 



Gt = 



fdr\F + (r)\* 
J drp 2 (r) 

and longitudinal magnetization, 

fdr\F z (r)\ 2 



G L = 



J drp 2 (r) 



(40) 



(41) 



for Q 



2ir x 237 Hz and 5 q = 0.2. 



2^90(90 + 2cin) 
This frequency is resonant with 2E < ^2 Q /h, and the fc ~ 
modes are expected to be excited according to the effec- 
tive Hamiltonian l|27p. From Fig. [1] (a), we find that the 
time constant for the exponential growth of Gt is ~ 22 
ms. 

Solving the Heisenberg equation for the effective 
Hamiltonian with k = 0, 



fl2 



cir 



-i£o (b+,ob-,o ~ &+, &L, 0) 



we obtain 



b± (t) = S± e- tE o h/n cosh£ i + ^ Q e iE ° h t h sinh£ t, 

(43) 



(a) 




200 400 600 800 1000 
t [ms] 

(b) |F + |/n argF + Fz / n 

t [ms] ^^HHIIH 

"tei 







285 



319 



800 




FIG. 1: (Color) (a) Time evolution of the average squared 
transverse magnetization Gt (red curve) and longitudinal 



magnetization Gl (blue curve) 
10 7 , 



for n 



2.8 x 10 1 



N = 10', fi = 2^/q (q + 2an) ~ 2tt x 237 Hz, and <5 3 = 0.2. 
The inset magnifies the dashed square, (b) Profiles of |_F+|/n, 
argF+, and F z /n at t = 0, 285, 319, and 800 ms. The field 
of view is 100 fjxa x 100 (im. The red circles indicate the 
topological defects. 



where £0 = qo5 q \uoVo\V . Substitution of this solution 



(42) into Eq. (HD with 



9±1 



(r,t) =J2 [Mr)b±.k(t) 



(44) 
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gives 



F+(t) 



2n 2^ e 



k 

2n(\u \ + \vq\) 



(|«fc| + |«fc|) blJt) + b_^ k (t) 



(45) 



where in the second line we take only the k — compo- 
nent because of the exponential factor e^ '. The average 
squared transverse magnetization Gt thus grows expo- 
nentially with a time constant l/(2£o) — 20 ms, which is 
in good agreement with the numerical results of 22 ms. 

The spatial profile of the magnetization is shown in 
Fig- HI (b). The transverse magnetization first grows 
with long wavelength [the second row of Fig. [1] (b)], in 
which we can see a few topological defects (red circles). 
These spin vortices are generated through the Kibble- 
Zurek mechanism (33l . [36j . The transverse magnetization 
then exhibits an interesting concentric pattern [the third 
row of Fig. [I] (b)], which may be due to the nonlinearity. 
The magnetization finally breaks into complicated frag- 
ments and the longitudinal magnetization also begins to 
grow [the fourth row of Fig. Q] (b)]. 

We now discuss symmetry breaking in the DCE. The 
original Hamiltonian @ commutes with J F z dr and the 
system has spin- rotation symmetry around the z axis. 
In reality, however, local magnetization in the x-y direc- 
tions occurs in a ferromagnetic BEC breaking the spin- 
rotation symmetry spontaneously [34j . and we expect 
that the symmetry breaking also occurs in the present 
system. This is because the exact quantum state with 
spin-rotation symmetry is the macroscopic superposition 
of magnetized states in the x-y directions, which is there- 
fore fragile against external perturbations. Such symme- 
try breaking phenomena should also occur in the DCE 
of the electromagnetic field. For example, in the ID 
moving mirror problem Q, there is rotation symmetry 
around the axis of the ID cavity and the generated pho- 
ton state is a superposition of all polarizations, giving 
(E) = (B) — 0. When we measure the local electromag- 
netic field, however, we obtain a nonzero value as a result 
of the symmetry breaking. 

Figure [5] shows the magnitude of the Fourier transform 



of F. 



+ - 



F+(k) 



drF+{r)- 



(46) 



Initially only the k ~ modes grow as shown in Fig. [2] 
(a). Then, the momentum distribution becomes broad 
as shown in Fig. [2] (b). Specific wave numbers can be 
selectively excited using larger O. Figure [2] (c) shows 
the result for SI = 2tt x 358 Hz, which is resonant with 
2E^ /% for k = 0.5 jivar 1 . We can see the ring at this 
wave number. 

From the form of the effective Hamiltonian in Eq. (f2"Tf , 
generated magnons are expected to have quantum corre- 



a o 



□ 




o 


1 0,1 





FIG. 2: (Color) Magnitude of the Fourier transform of F+ at 
(a) t — 285 ms and (b) t = 800 ms. The parameters are the 
same as in Fig. [1] (c) Magnitude of the Fourier transform of 
F+ at t = 500 ms for Q — 2n x 358 Hz. The other parameters 
are the same as in Fig. [1] 



lations. Defining new operators, 



£>+,fc — —^[b+^k + b-^k 



B-,k — —^[b+^k — b-^k 



(47) 
(48) 



we can rewrite Eq. (|27|) as 



,k 



B 



12 



.k 



(49) 



where = qoS q \ukVk\V. We can clearly see that this 
Hamiltonian generates the squeezed state. 

Figure[3](a) plots the values of -B+ : o = (&+,o + &- : o)/v2 
at t = and 30 ms. The distribution of B + corresponds 
to the quantum fluctuation in Eq. (|47[) with fc = 0, which 
is expected to be squeezed. The values of b±,o are ob- 
tained from b± k = W[\v,k\ip±i(k) — |«fc|V>^i(fe)] with 
ijj±i(k) being the Fourier transform of i^±i(r). Each 
point in Fig. [3] (a) corresponds to a single simulation run 
and 1000 simulations are performed. Figure [3] (b) shows 
the variance of Re(B + fle~ l9 ), 

V(0) - {[Re(B +fi e- i6 )] 2 ) ave - (Re(B + , e- ie ))l vg . (50) 

For the initial state, the variance is isotropic, V{ff) = 1/4, 
from Eq. (f39|) . The small deviation of the solid curve from 
1/4 in Fig. [3] (b) is the statistical error. At t = 30 ms, 
the distribution of B + $ is clearly squeezed, and the max- 
imum and minimum values of V(&) are 1.09 and 0.057. 
For the present parameters, e 2 ^' ~ 4.38 at t = 30 ms, 
which is in good agreement with 1.09/0.25 ~ 0.25/0.057. 
Since 1.09 x 0.057 ~ 1/4 2 , the squeezed state is almost 
the minimum uncertainty state. 

In this section, we assumed the use of 23 Na atoms. 
Similar results can be obtained also for F = 1 87 Rb atoms 
with a ferromagnetic interaction (ci < 0), where q(t) 
must always be larger than |ci | \I7§ in order to suppress 
the spontaneous magnetization. 
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FIG. 3: (a) Distributions of B+, = (6+,o + &-,o)/v / 2 at t = 
and 30 ms obtained by 1000 simulation runs with different 
initial states produced by random numbers. The parameters 
are the same as in Fig. [T] (b) Variance V(9) of Re(B+ i oe _ie ) 
for the data in (a). The dotted line indicates V(ff) = 1/4. 



V. DISCUSSION AND CONCLUSIONS 

We now discuss the possibility of experimental obser- 
vation of the proposed phenomena. The magnetization 
profile F(r) as shown in Fig. Q] can be measured by 
spin-sensitive phase-contrast imaging in a nondestructive 
manner [IS 0|, from which Gt, Gl, and F(k) are ob- 
tained. The squeezing of the field as shown in Fig. [3] can 
also be observed. From Eq. (|4"5|) , we can measure both 
b*^ k + b-.-k and b*^_ k — 6-,-fe, and therefore we obtain 
B±,fe. In order to assure that the observed magnetization 
is definitely due to amplification of the vacuum fluctua- 



tion, we must prepare the appropriate initial state. In 
principle, the initial magnon vacuum state can be pre- 
pared as follows. First, a BEC is prepared in the m = 
state at a sufficiently strong magnetic field (q ^> ci'Fq), 
in which the m = state is almost the ground state. The 
residual atoms in the m = ±1 states must be eliminated 
completely. Then, the magnetic field is adiabatically de- 
creased to the desired strength, which gives a magnon 
vacuum state satisfying Eq. (fl"0|) . 

In conclusion, we have studied the magnon excitation 
in a spinor BEC by a driven external magnetic field and 
have demonstrated the close analogy of this phenomenon 
with the DCE. The present system is suitable for study- 
ing the DCE of quasiparticles, since the time-dependent 
magnetic field applied to the Bogoliubov ground state 
amplifies only the vacuum fluctuation, keeping the "clas- 
sical fields" constant in the Bogoliubov approximation. 

We numerically demonstrated magnon excitation in a 
spinor BEC using the mean-field GP equation, in which 
the vacuum fluctuation is taken into account by the initial 
random noise. We have shown that the oscillating exter- 
nal magnetic field resonantly amplifies the vacuum fluc- 
tuation, leading to magnetization of the system (Fig.[T]). 
The Fourier transform of the excited field reveals that 
the specific wave number k can be selectively amplified 
(Fig. [2]). The excited quantum field is squeezed (Fig. [3]) 
as in the DCE of photons. The growth of magnetization 
in Fig. [T] (a) and the degree of squeezing in Fig. [3] can be 
well described by the effective Hamiltonian in Eq. (|2"T)) . 

A spinor BEC is thus a good testing ground for the 
DCE and our proposal is feasible with current experi- 
mental techniques. Study of the DCE of quasiparticles 
may serve as a stepping stone to the observation of the 
DCE of photons. 
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